Size measurement of dynamically generated resonances with finite boxes 
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The structure of dynamically generated states is studied from a viewpoint of the finite volume 
effect. We establish the relation between the spatial size of a stable bound state and the finite 
volume mass shift. In a single-channel scattering model, this relation is shown to be valid for a 
bound state dominated by the two-body molecule component. We generalize this method to the 
case of a quasi-bound state with finite width in coupled-channel scattering. We define the real- 
valued size of the resonance in a given closed channel using the response to the finite volume effect 
on the channel. Applying this method to physical resonances we find that A(1405) and /o(980) are 
dominated by the KN and KK scattering states, respectively, and that the distance between KN 
(KK) inside A(1405) [/o(980)] is 1.8-1.9 fm (2.7-2.9 fm). 

PACS numbers: 11.80.Gw, 14.20.-c, 14.40.-n, 11.30.Rd 



I. INTRODUCTION 

There are several hadrons which are expected to have 
some exotic structures (exotic hadrons), and clarifying 
structures of these exotic hadrons is one of the impor- 
tant tasks for the study of the strong interactions 
A classic example of exotic hadrons is the hyperon res- 
onance A(1405), which is the lightest baryon with spin- 
parity J p = l/2~ although containing one strange quark. 
This resonance has been considered as a quasi-bound 
state of the KN system 0, Q owing to the strongly 
attractive KN interaction in the 1 = channel. An- 
other example is found in the lightest scalar meson nonet 
[/o(600) = a, k(800), /o(980), and a (980)], which ex- 
hibits an inverted spectrum from the naive expectation 
with the qq assignment. There are several attempts to 
explain this anomaly, e. 3., multiquark configurations for 
the scalar nonet Ij, [5| and KK molecules for Jo (980) 
and a (980) @, [J Recently A(1405) and the light- 
est scalar mesons are successfully described by coupled- 
channel chiral d yna mics (chiral unitary approach) in 
meson-baryon [8rfl2| and meson-meson |la - tl5| scatter- 
ings, respectively. 

One of the characteristic features of exotic hadrons 
is the spatial size, because one expects larger size of 
hadronic molecules than ordinary hadrons. However, in 
general, candidates of exotic hadrons are not in ground 
states but resonances with finite decay width. Because 
of the decay process, mean squared radius of a resonance 
is obtained as a complex number whose interpretation is 
not straightforward [ljl [TtJ ■ To overcome this difficulty, 
we recall the finite volume effect on bound states. It has 
been shown in Refs. [HI, [l9| that the binding energy in- 
creases when a bound state of two particles is confined 
in a finite box with periodic boundary condition. The 



reason is that the wave function of the bound state in 
the box penetrates to the adjacent box and hence the ex- 
pectation value of the potential energy grows negatively 
This means that the finite volume effect is closely related 
with the spatial structure of the bound state. 

Motivated by these observations, in this study we aim 
at establishing the relation between the finite volume ef- 
fect and the spatial size of both stable bound states and 
unstable resonance states. Firstly, we consider a stable 
bound state in single-channel scattering where the size 
of the state is well defined. We develop a method to 
evaluate the size from the finite volume effect, and ex- 
amine its validity using a dynamical scattering model. 
This method is straightforwardly generalized to a bound 
state in coupled-channel scattering. In this case, the size 
of the bound state is defined for each channel, which can 
be estimated by the finite volume effect on the channel of 
interest, with the other channels being unchanged. Next 
we extend this method to a resonance state in coupled- 
channel scattering, and estimate the size of the state in 
closed channels. As applications to physical states, we 
examine the coupled-channel models for A(1405) and the 
scalar mesons a, /o(980), and a (980) to elucidate their 
structures. 

This paper is organized as follows. In Sec. [TT] we for- 
mulate the size measurement of (quasi-)bound states us- 
ing the finite volume effect, and introduce a dynamical 
scattering model. In Sec. IIIII we examine the validity of 
our strategy using the finite volume effect in the case 
of single-channel bound state, and apply the method to 
physical hadron resonances. Section [IV] is devoted to the 
conclusion of this study. 



II. FORMULATION 



A. Size measurement with finite volume effect 
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Here we present the basic idea to determine the spa- 
tial size of the (quasi-)bound state from the mass shift 
due to the finite volume effect. Let us first consider 
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FIG. 1: Schematic figure of (quasi-)bound state poles and 
continuum states in the complex plane of energy W. (a) A 
bound state in single-channel scattering, (b) a bound state 
in coupled-channel scattering, and (c) a quasi-bound state in 
coupled-channel scattering. 



the simplest bound state with mass Mb in the 

single-channel scattering of particles with masses to and 
to' < to [see Fig. [T](a)]. In nonrclativistic quantum me- 
chanics, the mean squared distance of two particles in an 
weakly bound state can be read off from the tail of the 
wave function as [l7| 



(r 2 ) NR 



4/iB E ' 



(1) 



with the reduced mass \x = mm' /{m + m') and the bind- 
ing energy = m+m' — Mb. At first glance, the spatial 
size of the bound state seems to be solely determined by 
the mass of the bound state. However, it is implicitly 
assumed in Eq. ([!} that the bound state is completely 
described by the model space of the two-body scattering. 
In field theory, in addition to the scattering state, there 
can be a "bare state" contribution (elementarity) whose 
fraction in the physical bound state is expressed by the 
wave function rcnormalization constant Z 12(1 I21L 1 The 



1 Strictly speaking, Z represents the contribution from those other 
than the present model space. Although this is not necessarily 
an elementary particle, for simplicity, we call Z the "bare state" 
contribution in this paper. Also we consider a pointlike bare 
state which does not contribute to the size of the bound state. 
This is in line with Ref. [lgl ] where the finite volume effect is 
attributed to the modification of the loop of scattering states. 



size of the bound state, which stems from the scattering 
state contribution, should then be given by subtracting 
the bare state contribution as 



1 - Z 



(2) 



where the factor 1 — Z is called the compositeness. As 
shown in Refs. [20rl22| . the compositeness is related with 
the coupling constant of the physical bound state to the 
two-body scattering state g: 



l-Z = -g 2 G'(M 2 ), G'(s) = 



dG(s) 
ds 



(3) 



where G is the two-body loop integral as a function of 
squared energy s = W 2 to be specified below. Thus, the 
size of the bound state is expressed in terms of the mass 
Mb and its coupling to the scattering state g. 

At this point we make use of the finite volume effect 
on the mass of a bound state studied in Refs. [TH, [l9j . 
When a bound state is put in a periodic finite box of 
size L, the mass Mb is shifted to M-q(L), and the mass 
shift AA/b(£) = AIb(L) — Mb is related to the coupling 
constant g. As shown in Appendix [XJ the leading con- 
tribution to the mass shift formula in the present case is 
given by 



AM B {L) = - 



3g 2 



8irM 2 L 



exp[-^L]+C9(e-^ i ), 



/' 



y/-A(JV/g, to 2 , m' 2 ) 
2M^ 1 



(4) 
(5) 



where \(x, y, z) = x + y + z — 2xy — 2yz — 2zx is the 
Kallen function. An important point is that the leading 
contribution depends on g 2 and Mb in the infinite vol- 
ume. Namely, we can read off the physical coupling con- 
stant g from the L dependence of the mass of the bound 
state. Equations @ and Q show that the bound state 
has large size when the coupling g is large. This fact 
is intuitively understood in Eq. Q by the g 2 factor; a 
bound state with large size (and thus the large coupling 
g) has strong finite volume effect. In other words, the 
structure of the bound state is quantitatively reflected in 
the finite volume effect. 

Equation (@| provides us an alternative strategy to es- 
timate the size of a bound state. Suppose that we are 
able to calculate the mass shift AMb(L) of the bound 
state in finite volume at large L. In this case, using 
the mass shift formula (UJ, we determine the coupling 
strength gpv- The mean squared distance is then evalu- 
ated as 



/FV = 



-9lyG'{M 2 ) 



(G) 



If the L dependence of the mass shift is correctly fitted 
by Eq. (j4j, we expect (?fv — > g and (r 2 )py — > (r 2 ). In 
this way, the size of the bound state is related with the 
finite volume mass shift at large L. 
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The virtue of this new approach will become clear when 
the argument is extended to the quasi-bound state with 
finite width. To this end, we begin with the case of 
a bound state in multi-channel scattering as shown in 
Fig. [T] (b). Labeling the scattering states by the suffix i, 
we can decompose the bound state wavefunction into the 
bare state contribution (Z) and the contribution from 
the scattering state in channel i (JQ ) [22[ , which are nor- 
malized as 



1 =Z 



X; 



i 



(7) 

(8) 



In this case, the mean squared distance of the bound 
state in channel i is defined as 



(9) 



with gi being the coupling constant to channel i, B^i = 
m, + m'i — Mb and fii = mi'm'J (m^ + m^). Conceptually, 
(r 2 ),; corresponds to the size measured by the probe cur- 
rent which exclusively couples to the component in chan- 
nel i. Then the coupling constant (^fv can be extracted 
from the mass shift by putting only the channel i in the 
finite box with size L and keeping the other channels un- 
changed. Substituting this coupling constant .g^FV into 
Eq. ((9]), we obtain (r 2 )j,FV- 

We can further extend this argument to a quasi-bound 
state with finite width. Consider a system with two cou- 
pled channels in which the higher energy channel has 
a bound state when the transition potential is switched 
off. The bound state acquires a decay width through the 
channel coupling to the lower energy channel, which is 
called a quasi-bound state or a Fcshbach resonance [see 
Fig- HI (c)]. In this case, the pole of the resonance locates 
in the complex energy plane at total energy W = W po ie- 
If the channel coupling is not strong, the imaginary part 
of the pole position is small and we can identify the real 
part as the "mass" of the state, Mb = Re[W po ic]- Ap- 
plying the same procedure, we determine the coupling 
constant g^FV from the L dependence of the real part 
of the resonance pole when the channel i is put in the 
box. Substituting it in Eq © , we estimate the size of the 
quasi-bound state. The binding energy as well as the loop 
function are evaluated at this energy Mb = Re[W po i e ]- 
Note that this is only applicable to the closed channels, 
namely, the resonance pole should be located below the 
threshold of channel i. If we put an open channel in 
the finite box, the continuum state of that channel is dis- 
cretized and we cannot perform the analytic continuation 
to the complex energy plane. 

It is important that this procedure gives a real-valued 
(r 2 )i of the quasi-bound state, since the coupling ex- 
tracted from Eq. (|4]) is a real number. In general, it 
is known that the size [l(| 53] and compositeness [22| 
become complex in the case of resonances, which are 
difficult to interpret. The strategy presented here can 



provide an alternative way to investigate the structure of 
resonances. In the following, we introduce a hadron scat- 
tering model together with finite volume effect, in order 
to examine the size measurement with finite box. 



B. Coupled-channel scattering model and finite 
volume effect 



Here we formulate a model to describe stable bound 
states and unstable resonance states along the line with 
Ref. [l(J. We prepare a coupled-channel interaction ker- 
nel Vij and evaluate the scattering amplitude T,j by the 
Bcthe-Salpctcr equation in its factorized form: 



Tij(s) 



s2v ik G k T kj = J2(l-VG)-\ k V kj , 



(10) 

where indices i, j, and k represent the scattering chan- 
nels, s = W 2 is the squared center-of-mass energy of the 
scattering system. The explicit form of will be given 
in the next section. G k is the loop integral, 



Gfc(s) 



i 



i 



(2tt) 4 q 2 - m 2 + it (P - q) 2 - m' 2 + ie 

(11) 

with m k , m' k , and = (W,0) being the masses of the 
particles in channel k and the four-momentum of the two- 
body system, respectively. Using the dimensional regu- 
larization, one can rewrite the loop integral as 



G fc (s) 



1 



16tt 



afe(Mrog) + m-2— 

Mreg 




mi 



(12) 



with the rcgularization scale /U reg , the subtraction con- 
stant a k , and X k = A(s, m 2 ., m' 2 ). We note that the 
rcgularization scale and the subtraction constant are not 
independent, and the subtraction constant is a single pa- 
rameter of the loop function in each channel. 

Bound states and resonance states appear as poles in 
the scattering amplitude as 



Tij (s) 



9i9j 



Spolc 



(13) 



where the background term T^ G is chosen to make the 
product gigj energy independent. The constant gi can 
be interpreted as the coupling strength of the state to 
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the channel i. The pole position s po i c is a solution of the 
equation, 



det(l - VG) = 0, 



(14) 



which is simplified as V^ 1 = G in the single-channel case. 
A stable bound state is represented by a pole on the real 
axis of the first Riemann sheet below the threshold, while 
an unstable resonance state corresponds to a pole in the 
complex lower-half plane of the second Riemann sheet 
above the threshold. 

In this model, we have a relation, 



E 



dGi 



® as 



1, (15) 



which is the generalized Ward identity [17J, 12J, |23J . With 
Eq. ©, we can identify the first term as the sum of the 
contributions from hadronic composite states. It follows 
from Eq. (J7J that the bare pole contribution Z is ex- 
pressed by the second term as 



dV, 



(16) 



In Rcf. [22j Z is shown to be exactly the bare pole con- 
tribution for stable bound states, and thus for unstable 
resonance states the system is expected to have less com- 
positeness as Z approaches unity. 

Next we consider the finite volume effect in a spa- 
tial box with size L. In general, the finite volume effect 
appears as the discretized momentum in the loop func- 
tion [24|]. The finite volume effect in the present model 
has been discussed in Refs. [25|, [26| by using discretized 
momentum loop function Gi, 



dq° 



1 



2tt 



(P - 1? 



2-Kn 



(17) 



instead of the loop integral Gi in Eq. (fTQj) . Here, we 
evaluate Gi with the dimensional regularization following 
Ref. p6| by extracting three-dimensional integral from 
Gi and replace it with the summation with discretized 
momentum, which results in, 

Gi(s) -> Gi(s) = Gi(s) 



3 Ii(s,q) 



with 



h(s,q) 



1 UJj + u)[ 

2oJiOj' i s — (cjj + ui^) 2 ' 



(18) 



(19) 



r 



(20) 



It is known that, with finite cut-off q ma , x , Gi exhibits 
oscillations which gradually vanish as g ma x goes to in- 
finity [26]. This oscillation is caused by the summation 
over the discretized momentum, which is not a continu- 
ous function of q max - The absolute value of the integrand 
7i(s, q) decreases for large q, so the discontinuity becomes 
small with large g max . In order to make convergence with 
respect to the oscillation, following Ref. [26| , we will take 
averaged value of Gi within range q max £ [2 GeV, 4 GeV] 
in the numerical calculation. 

When we put all the channels in the finite box, mo- 
menta of the scattering states above the threshold are 
also discretized, and the eigenenergies are constrained by 
the pole condition [25| 



det(l - VG) = 0, 



(21) 



which again reduces to V~ x — G in the single-channel 
case. Note that for bound state poles below the thresh- 
old, Eq. (|T4)) and Eq. ((21]) are the same condition with 
different loop function. For the application to the quasi- 
bound state, we will use the loop function with channel 
i in the finite box as 



(G 1 



G = 



\ 



G-, 



G, 



\ 



(22) 



In this case, if the energy is smaller than the threshold 
of channel i, W < rrii + m£, the scattering amplitude 
T = (V^ 1 — G) _1 is a continuous function of W, and the 
resonance pole can be searched for through the analytic 
continuation of the amplitude in a usual manner. 



III. RESULTS 
A. Size of bound states in single-channel scattering 

Now let us consider stable bound states in single- 
channel scattering and see how they behave in the finite 
volume. In Sec. Ill Al we have presented two methods to 
calculate the size of the bound state. The mean squared 
distance in Eq. © is obtained from the residue of the 
pole, and that in Eq. © is evaluated by the finite vol- 
ume effect. In addition, the corresponding mean squared 
distance can also be calculated by using the response to 
an external probe current as shown in Ref. [ItJ • Compar- 
ing the results from different methods, we examine the 
validity of the size estimation. 

In this subsection, with the KN system in mind, 
we choose the masses in the scattering state as m = 
938.92 McV and m' = 495.67 MeV, respectively. We use 



TABLE I: Properties of bound states in cases I and II. The 
results in column "Pole" are calculated from the residue of the 
pole in the models, and those in "Finite Volume" are obtained 
by using the finite volume effect. 



Case I 


Pole 


Finite Volume 


9 
Z 
1 - Z 


5.415 GeV 

1 


5.170-5.772 GeV 
-0.136-0.089 
0.911-1.136 




Case II 


Pole 


Finite Volume 


9 
Z 
I - Z 


1.383 GeV 
0.935 
0.065 


1.679-1.994 GeV 
0.864-0.904 
0.096-0.136 



the natural subtraction constant a = —1.95 with the reg- 
ularization scale /x reg = 630 MeV, which is obtained to 
exclude explicit pole contributions from the loop inte- 
gral [l?]] ■ For the interaction kernel V we consider two 
types. One is the constant interaction, 

Vi = v , (23) 

with the energy independent parameter vq (case I). The 
other interaction consists of a bare pole term, 

V U = (24) 

s - s 

with two parameters go and sq which are constrained by 
9o = f° r simplicity (case II). The parameters vq and 
.9o = y/so are nx °d so as to produce a bound state with 
binding energy Be = 10 MeV in both cases, and as a 
result we have v = 124.9 and go = y/so = 1-430 GeV. 

Properties of the bound states in two cases I and II 
are summarized in the second column of Table HI The 
coupling constant g is calculated from the residue of the 
bound state pole as in Eq. (|13p . and the bare pole con- 
tribution Z is obtained by Eq. (|16p. As one can see 
from Table HI the bound state by the constant interac- 
tion in case I has Z = [2a |, which can be understood 
by Eq. (|16p. On the other hand, the bare pole potential 
creates large elementarity Z = 0.935 in case II. Purely 
elementary state with Z = 1 can be obtained by taking 
the limit go — > and so — > M| with fixed Be- 

With the obtained compositcness 1 — Z and Eq. (J3J, 
we calculate the size of the bound state as shown in 
the third column of Table [TTJ Here we also calculate the 
mean squared distance using probe method developed 
in Ref. [171 ]. in which the external probe current is cou- 
pled to the particles in the scattering state and the mean 
squared distance is obtained from the form factor. The 
results are shown in the second column of Table [TT1 Com- 
paring two cases, we observe that the bound state in 
case I has large size of \J (r 2 ) ~ 1.9 fm compared to the 
hadronic scale < 1 fm, whereas in case II the size of the 
bound state is \J (r 2 ) ^0.5 fm. This is because only the 
two-particle cloud can contribute to the mean squared 
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TABLE II: Sizes of bound states in cases I and II with several 
methods. 







Probe [17.1 


Eq. © 


Eq. © 


Case I 


<r 2 > [fm 2 ] 


3.474 


3.001 


2.735-3.410 


Case II 


0.227 


0.195 


0.288-0.407 





-0.002 
-0.004 



G 

G(L = 4.0fm) 


1 
f 
* 
t 
1 


1 ■ 
i 

\ . 


G(L = 5.()fm) 

G(L = 6.()fm) 


t 
* 
1 
1 
t 


\ 

\ 




* 
* 
* 




V5 1 


* 
f 
* 

* 
* 
* 
* 






/ 






\/\\ 
/•- . \ \ 





1400 1405 1410 1415 1420 1425 1430 1435 1440 

W [MeV] 



FIG. 2: The loop integrals in infinite and finite volume, G 
and G, and inverse of the interaction Vi and Vn as functions 
of W = y/s. Vertical dotted line represents the threshold of 
the two-body system, W — m + m! . 



distance. These results indicate that the spatial size of 
the bound state is not exclusively determined by its bind- 
ing energy, and the magnitude of the coupling constant 
is closely related with the internal structure, as discussed 
in Sec. UTAl 

For later convenience, the loop integral G and the in- 
verse of the interaction kernels V^ 1 and V^ 1 are plotted 
as functions of W = y/s in Fig. [21 Below the threshold, 
the loop function G is real. In this figure, the intersec- 
tion point of G and V~ l corresponds to the mass of the 
bound state Mb according to Eq. (fl"4")) . In both cases the 
intersection appears at Mb = 1424.6 MeV with the bind- 
ing energy Be = 10 MeV with the adopted parameters. 
An important point to note here is that the energy de- 
pendence of two interaction kernels is very different from 
each other. While V^ _1 is completely flat, V^ 1 is almost 
vertical with steep slope. In the limit of go — > and 
So — > Mg, the slope becomes completely vertical. This 
difference of the interaction kernel will be crucial to the 
finite volume effect on the bound states. 

Then let us take into account the finite volume effect 
by replacing the loop integral G with that in finite vol- 
ume G. Behavior of G is also plotted in Fig. [5] with box 
sizes L — 4.0, 5.0 and 6.0 fm. Because of the pole con- 
dition (|2"Tj) , the mass of the bound state in finite volume 
Mb(L) corresponds to the intersection point of G and 
V^ 1 . From Fig. [5[ one observes that in both cases the 
mass of the bound state decreases when the box size L 
decreases. However, L dependence of the mass of the 
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the factor — 3exp(— ^,L)/(8ttM^L) as, 



bound state is quantitatively different in two cases. The 
fiat (steep) energy dependence of Vf 1 {V{[ ) results in 
the strong (mild) L dependence of the bound state mass 
in finite volume. Different L dependence of the mass shift 
in two cases is understood by this geometric argument. 

To compare with the mass shift formula we plot 
in Fig. E] the mass shift AM B (L) = M B (L) - M B as a 
function of L. From this figure, we observe the decrease 
of the mass for the smaller box size L in both cases I 
and II. Furthermore, one can see the rapid decrease of 
the mass in case I compared to that in case II. This can 
be interpreted as the consequence of the loose binding of 
the system (large mean squared distance) in case I. 

Using the coupling constant obtained from the pole 
residue, we can predict the mass shift AM B by Eq. (|4]) 
which is plotted in Fig. [3] With large L, the formula 
well reproduces the mass shift, but some deviation be- 
comes evident in smaller L region, especially for the case 
II. This means that higher order corrections on the mass 
shift formula is necessary to describe finite volume effect 
of the bound state . In fact, since the coupling g is small 
in the case II, it is reasonable that the higher order cor- 
rection to the mass shift formula is more important than 
the case I. 

Let us extract the bound state properties by using the 
finite volume effect in the procedure of Sec. Ill Al Fitting 
the mass shift by the formula (|4|, we evaluate the cou- 
pling strength <?FV; the bare pole contribution Z-pw and 
mean squared distance (r 2 )FV- In this study we take the 
following strategy to evaluate the coupling strength gpy. 
Namely, since the mass of the bound state is expected 
to change according to Eq. (0| at the leading order, the 
coupling strength gpy can be extracted by the fraction 
of the mass shift obtained in our model, AM B (L), and 



5fv 



i 



AMb(I) 



87rAf|L 



exp [— fiL] 



(25) 



with 



2Mv 



(26) 



This gpy depends on the box size L especially in small 
L region where the higher order contributions are not 
negligible. Nevertheless, we expect that gpy m Eq. (f25|) 
becomes almost flat in the region where the mass shift is 
dominated by the leading order contribution. In Fig. |4] 
we plot i?fv in Eq. (|25[) as a function of the box size L for 
both cases I and II. From the figure, we can see that <?fv 
in case I is fairly flat at L ~ 3 fm, while it rapidly changes 
below ~ 2 fm due to the higher order contributions to the 
mass shift. On the other hand, in case II <7fv increases 
without flat regions as L decreases down to ~ 1 fm. This 
reflects that the higher order contributions to the mass 
shift are not negligible even at L ~ 4 fm. Here in order to 
suppress the higher order contributions we take fit ranges 
[2 fm, 8 fm] and [4 fm, 8 fm] in case I and II, respectively, 
and adopt gpy in these ranges as the coupling strength 
from the finite volume effect. The adopted values of gpy 
is shown as bands in Fig. [4] for both cases I and II. The 
results are summarized in Tables HI and HT1 As a result, we 
qualitatively reproduce the structure of the bound state. 
Especially the properties of the bound state in case I 
are reproduced within ~ 10% accuracy. This indicates 
that the size measurement with the finite volume effect 
is a powerful tool to clarify the structure of bound states 
which have dynamical origin. 
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B. Application to physical resonances 

In the previous subsection we have developed a method 
to estimate the size of the bound state by using the finite 
volume effect. One of the important features of our pro- 
cedure is the applicability to Feshbach resonance states 
with finite widths as discussed in Sec. Ill Al Furthermore, 
one can obtain real-valued size of the resonance states 
with respect to a closed channel. In this subsection we 
use this method to discuss the structure of physical res- 
onance states from the finite volume. 

Let us discuss A(1405) in KN-ir'E-rjA-K'E. coupled- 
channels and a, /o(980), and ao(980) scalar mesons in 
■RTT-'K-q-KK coupled-channels, assuming the isospin sym- 
metry. These resonances have been studied in chiral 
unitary approach [||-[l5j], which is now elaborated using 
next-to-leading order chiral interactions with recent ex- 
perimental data (28l - l30j . To concentrate on the size es- 
timation with the finite volume effect, here wc utilize 
simplified models with leading order interactions as fol- 
lows. For A(1405) we employ the Weinberg- Tomozawa 
term as the interaction kernel, 

Mi-Mfiy/tMiMj, (27) 



with Mi being the baryon mass in channel i, f the meson 
decay constant, and Cy the Clebsch-Gordan coefficient 
which is determined by the SU(3) group structure of the 
interaction, 



C ; 




3/\/2 \ 
4 

-3/\/2 

3/2 -3/V2 3 J 



(28) 



where i = 1, 2, 3, and 4 denote the KN, irY,, r/A, and 
ifS channels, respectively. The meson decay constant 
is / = 1.123/^ with f n = 93.0 MeV. The subtraction 
constant is a\ = —1.84, = —2.00, 03 = —2.25, and 
G14 = —2.67 with the regularization scale ^ rcg = 630 MeV 
in all meson-baryon channels [3l|. For the scalar meson 
case, we take the lowest order s-wave meson-meson in- 
teraction in chiral perturbation theory as the interaction 
kernel, namely, 



V u = 



mj - 2s 
~2P ' 



Via = V21 = - 



\/3s 3s 

47T 22 = ~4/f 

(29) 

for the 1 = channel with i = 1 (2) for irir (KK) channel, 
and 



Vll = -7777, = V21 = 



3/ 



(9s— m^— 3m 2 n — 8m 2 K ), 



V22 = - 



4/r 



(30) 
(31) 



for the 1=1 channel with i = 1 (2) for irrj (KK) channel. 
Here we use the pion decay constant = 93.0 MeV. The 



TABLE 111: Properties of resonances in A" N-n^-i/A-KE scat- 
tering and nn-nrj-KK scattering. Pole positions (Wpoic), cou- 
pling constants (g), and decomposition into scattering states 
(X) and bare pole contributions (Z) are shown. 





A(1405), higher pole 


A(1405), lower pole 


TAT" 


14zD.z — ID. (1 MeV 


1 inn c: aft o„' A/T^w 
loyU.0 — OO. Zl MeV 


9kn 


6.021 + 2.157; GeV 


-2.689 + 4. 146; GeV 




-1.187 - 3.826; GeV 


6.255 - 4.096; GeV 




3.630 + 0.502; GeV 


0.030 + 1.980; GeV 


9KS 


-0.352 - 0.928; GeV 


1.220 - 1.184; GeV 


Xkn 


0.994 + 0.048; 


-0.213 - 0.125; 




-0.047 - 0.151; 


0.365 + 0.526; 


X V A 


0.052 + 0.012; 


-0.014 + 0.003i 




-0.002 + 0.002; 


0.000 - 0.007; 


z 


0.002 + 0.089; 


0.863 - 0.397; 



W pole 

9kk 

X n7 z 

z 



471.3 - 181. 0; MeV 
1.837 - 2.305i GeV 
0.802 - 1.161; GeV 
-0.163 + 0.347; 
-0.005 - 0.009; 
1.169 - 0.339; 



/o(980) 



Wpole 
Qtyty 

9kk 
X nn 

z 



987.1 - 17.7; MeV 
-0.483 + 1.428; GeV 
3.909 + 1.316; GeV 
0.014 + 0.010; 
0.739 - 0.107z 
0.246 + 0.097; 



a (980) 



Wpole 

9kk 

A71-77 

X-kk 
Z 



979.4 - 53.4; MeV 
-2.941 + 0.777; GeV 
4.579 + 0.481; GeV 
-0.059 + 0.104; 
0.377 - 0.288; 
0.682 + 0.184; 



subtraction constant is fixed at 



T with the regular- 



ization scale fi Teg = 1.325 GeV in all meson- meson chan- 
nels, which corresponds to the three-dimensional cut-off 
1.092 GeV [H|. 
With these interaction kernels, we obtain two reso- 
nance poles in the meson-baryon scattering amplitude 
below the KN threshold, both of which are associated 
with A(1405) [H, HII- I n meson-meson scattering, we 
find two poles in I = and one in I = 1 below the 
KK threshold, which are interpreted as er, Jo (980), and 
ao(980) mesons, respectively. Properties of dynamically 
generated resonances are summarized in Table IIIII The 
higher pole of A(1405) is expected to originate from the 
KN bound states (33[, and in fact the magnitude of the 
KN component X^- N is much larger than the others. In 
the scalar meson case, we have a meson with very large 
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-180 
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Re[W pole ] [MeV] 
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FIG. 5: Behavior of resonance pole positions for two A(1405) 
(a) and three scalar mesons (b) with several box sizes L. Here 
filled symbols indicate pole positions in infinite volume and 
open points are plotted in interval 0.5 fm with respect to the 
box size L from L — 3.0 fm to L = 7.0 fm. Note that a pole 
is quite stable with respect to the finite volume effect on the 
KK channel and the pole for ao(980) disappears for box sizes 
smaller than 4.5 fm. 



width ~ 400 MeV together with / (980) and a (980)^. 
In the present setup, /o(980) is dominated by the KK 
component whereas do (980) shows large bare pole con- 
tribution Z. We note that the pole positions of /o(980) 
and ao (980) depend on the cut-off for the loop integral 
and with smaller /x reg they move above the KK thresh- 
old [l4|. In this case, the quasi-bound state picture for 
Jo (980) and ao(980) becomes unclear. 

Then let us take into account the finite volume effect. 
Since they are the closed channels for all the poles con- 
sidered here, we put KN and KK channels into finite 
boxes with the periodic boundary condition with other 
channels being unchanged. Behavior of the resonance 
pole positions with respect to the box size L is shown in 
Fig. [5J In the A(1405) case [Fig. [S] (a)], the higher pole 
moves to lower energies when the box size for the KN 




FIG. 6: Coupling strength evaluated by Eq. (]32|) as a function 
of L. The bands show the adopted values of <;fv in our study 
for higher pole of A(1405) and /o(980). 



channel is reduced. On the other hand, the lower pole 
stays around the original pole position even if the finite 
volume effect on the KN channel is taken into account. 
This indicates that the higher pole is largely affected by 
the modification of the KN loop and supports the sce- 
nario that this pole originates from the KN bound state. 
In the scalar meson sector [Fig. [5] (b)], a and ao(980) do 
not follow the expected mass shift formula; a is quite 
stable with respect to the finite volume effect on the KK 
channel and the shift of the pole position is less than 1 
MeV. The do (980) pole disappears for box sizes smaller 
than 4.5 fm. On the other hand, the pole position of 
/o(980) shows strong L dependence and moves to lower 
energies for smaller box size L. This implies large KK 
component inside /o(980), which is not prominent for a 
and a (980). 

We next estimate the size of dynamically generated 
resonances with the procedure developed in Sec. Ill Al In 
our approach, we identify the real part of the pole posi- 
tion as the mass of the state. This means that we expect 
a downward shift of the real part of the pole position 
in finite volume. However, this is valid only when the 
coupled-channel contribution to the pole under consider- 
ation is so small that we can neglect the mass shift from 
the coupled-channel effect. In fact, the poles for a and 
do, as well as the lower energy pole of A(1405), do not ex- 
hibit the downward mass shift in finite volume. We then 
conclude that these states are not dominated by the KN 
(KK) component, in agreement with the results in Ta- 
ble [Hi] Therefore, we here consider the properties of the 
higher pole of A(1405) and /o(980) resonance with re- 
spect to the KN and KK component , respectively. We 
first fit the coupling strength g^ N FV {qkk fv) to the L 
dependence of the real part of the pole position of the 
A(1405) [/o(980)], and then evaluate the mean squared 
distance between KN (KK) in A(1405) [/ (980)]. As 
in the case of the bound state, we extract the coupling 
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TABLE IV: Properties of A(1405) and /o(980) with finite vol- 
ume effect. 



A(1405), higher pole 


B E 


= m N + m K — Re[VK P oie] 


8.354 MeV 




9kn,fv 


4.859-5.152 GeV 




Xkn.fv 


0.890-1.001 




( r }kn,fv 


3.198-3.596 fm 2 




/o(980) 


Be 


= 2m K - Re[W po ie] 


4.218 MeV 




9kk,fv 


2.822-3.053 GeV 




X KK, FV 


0.783-0.916 




{ r }kk,fv 


7.287-8.526 fm 2 



strength g F y by, 



9kn(kk),fv 



Re[W p olc - Wpole 



\ 87r(Rc[IF pole ]) 2 L 



exp [-Hkn(kk) l ] 
(32) 



We can compare the present result with previous cal- 
culations of the KN distance in A(1405). In Ref. fjj\ . 
the mean distance of the KN bound system at 1424 MeV 
is calculated by the probe method, which leads to 

VTr^y =1.69 fm (probe method, B ~ 11 MeV). 

In Ref. [34[ , the mean distance is calculated by the effec- 
tive single-channel potential developed in Ref. [33j | . The 
results are 

^/(r 2 ) =1.72 - 1.99 fm (potential, B ~ 10 - 13 MeV). 

Comparing with these results, the present method gives 
consistent values of the root mean squared distance be- 
tween KN, 1.8-1.9 fm. However, we note that, since the 
mean squared distance in the present method is propor- 
tional to 1/Be and the Feshbach resonances considered 
here are loosely bound states with small binding energy 
< 10 MeV, the size of the state would be sensitive to 
the small shift of the binding energy of several MeV. In 
this way, we find that the precise determination of the 
pole position of A(1405) and /o(980) is important for the 
quantitative study of the structure of the resonance. 



^KN(KK) 



^-A((Re[M/ pole ]) 2 , 



K> 



l N(K) 



2Rc[^ po i, 



(33) 



for higher pole of A(1405) [/ (980)]. 2 Here W po ic and 
Wpoic are resonance pole position in the complex energy 
plane in infinite and finite volume, respectively. In Fig. [6] 
we plot g-pv in Eq. (|32j) as a function of the box size L. In 
the figure, we observe rapid change of <7fv also in large L 
region. This means that the pole of the resonance states 
does not simply move downward in large L region due 
to the coupled-channel effect. Nevertheless, we observe 
fairly fiat <7fv i n range [2 fm, 5 fm] for the higher pole of 
A(1405) and [3 fm, 5 fm] for / (980). Hence, we adopt 
gvy in these ranges as the coupling strength from the 
finite volume effect. The adopted values of $fv are shown 
as bands in Fig.[5]for both A(1405) and /o(980). 

The coupling constants and the estimated sizes are 
summarized in Tabic IIV1 As one can see from the ta- 
ble, these resonances are dominated by the KN (KK) 
component with large spatial extent. In addition, the 
magnitude of the KN (KK) component X^ N (Xkk) 1S 
in fair agreement with that obtained on the pole posi- 
tion presented in Table UTTl Root mean squared distances 
are y^j = 1.8-1.9 fm for A(1405) and 2.7-2.9 fm for 
/o(980), which are larger than the typical hadronic scale 
< 1 fm. In this way, the A(1405) and / (980) can be 
interpreted as loosely bound Feshbach resonances. 



IV. CONCLUSION 

In this paper the structure of dynamically generated 
hadrons has been discussed from the viewpoint of the 
finite volume effect. We have presented a method to ex- 
tract the properties of a bound state in single-channel 
scattering using the finite volume mass shift. Introduc- 
ing a dynamical scattering model, we have shown that the 
coupling strength, compositeness, and spatial size of the 
bound state in infinite volume can be reproduced with 
good accuracy from the mass shift of the bound state in 
finite volume. 

This technique have been extended to a quasi-bound 
state with finite width in coupled-channel scattering, pro- 
vided that the width is small. We can estimate the spa- 
tial size of the component in a closed channel from the 
movement of the pole position along with the finite vol- 
ume effect on this channel. For an application to physical 
resonances, we have considered A(1405), a, /o(980), and 
ao (980) described in chiral unitary approach for coupled- 
channel hadron scatterings. Applying the finite volume 
effect on the KN and KK channels, we have found that 
the poles for the higher A(1405) and /o(980) move down- 
ward in finite boxes. This result indicates that A(1405) 
and /o(980) respectively have large KN and KK compo- 
nents. Fitting to the mass shift formula, spatial distances 
of KN and KK have been evaluated as 1.8-1.9 fm and 
2.7-2.9 fm for higher A(1405) and / (980), respectively. 



2 Note that the masses of K and K are the same, but they are 
distinguishable by the strangeness quantum number. We use the 
formula l|4")l. which differs by factor 2 from the one for identical 
particles in Ref. [T^l . 
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Appendix A: Mass shift of bound states in finite 
boxes 

In this Appendix we derive the leading contribution 
to the mass shift formula for bound states in finite box 
of the size L, Eq. gj), following Refs. [H [H]. Here we 
consider a bound state with mass Mb coupled with a 
two-particle system with masses to and to' < to. In this 
appendix, we work in the Euclidean space. We consider 
the small binding region as 



to' 2 <A/b < Tn 



(Al) 



In a finite spatial volume, the momentum of the two- 
particle system is discretized as q(L) = 2irn/L with 
n G Z 3 and the mass of the bound state is shifted 
to be Mb(L) 7^ Mb- Expanding the self-energy in fi- 
nite volume £(P) around P° = 1Mb, the mass shift 
AM B (£) = M B (L) - M B is given by 



AM B (L)=- 



1 



2Mb 
(iM B , 



t(P) - £(P)] + O(AMl) 



(A2) 
(A3) 



where £(P) is the self-energy in the infinite volume. 
While several diagrams contribute to the self-energy [l!| , 
the leading effect to the mass shift stems from the dia- 
gram shown in Fig. [7] The momcntum-discrctized loop 
integral can be expanded in powers of e~ lLm ' q ,m G 1? 
with the help of the Poisson summation formula, and the 
leading contribution can be obtained as 



E(P)-£(P) = J -|^2 



x TG m [(l-S)P + q]G m ,(6P-q)T 
+ 0{e-^ L ), (A4) 



where T is the three-point vertex function, G m is the 
propagator with mass to, and 



t 1 =- 



2M~ r : 



with X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx. 

The momentum fraction S > is chosen to maximize 
the analytic region of Im qi as follows. Here we consider 
that the particles of Mg and to have a conserved charge 3 
so that we can trace the line which connect the external 
Mb and to propagators of the vertex function T. We 
then use the same argument with Ref. [l!| to assign the 
momenta of the internal lines in the vertex function T. 
The conditions to avoid singularity are found to be 



{l-5)P± l -q 



< to 2 < M|, 



Im { SP ± -q 



< (to') 



'■\2 



By choosing 



M| + (to') 2 - to 2 
2M% 



the maximum analytic region for Im q± is obtained as 

< (Im qi f < 4/r 2 , 

The poles of the propagators as functions of q\ are given 
by 



i^Qo V 2 + l ( 2 - 2S)M B qo, (A5) 

(A6) 



q 2 ± +fi 2 - i25M B qo, 



for G m [(l-5)P+q] and G m <(8P-q), respectively. Mod- 
ifying the integration contour properly, we obtain two 
terms from these poles with the rest contributions being 
higher order corrections after the q\ integration: 



with 

h = 6i 

I[ = Qi 



t{P) - S(P) =h 0(e'^ L ) 
dqod 2 q± cxp(iLqi) 



(A7) 



(2tt) 3 2 Ql 

dqod 2 q± cxp(iLqi) 
(2tt) 3 2^ 



TG ml (8P-q)T 

'<Jl=9l 

rG m [(l-<5)i > + g]r 



9l=9i 



3 In the application to A(1405), baryon number is conserved for 
A(1405) and JV. For /o in the KK scattering, we can apply the 
formula in the equal mass case m = m' by Liischer jla ] except 
for the symmetric factor 1/2, which coincides with Eq. I A10I I . 
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where we have used the rotational invariance. The re- 
maining propagators have a pole in the complex q Q plane 
at 



5o 



=0. 



(A8) 



The leading contribution to the mass shift formula comes 
from this pole. To obtain the saddle point expression, 
we shift the go integration path from the real axis to 
Im q = -i(l - 5)M B (Im q = iSM B ) in h (I[). The 
pole contribution from % is picked up by /{ term. The 
leading contribution is then given by 



I[=6 



d 2 qj_ cxp(-L v /^ 2 + qj_) g 2 
(2tt)3 

3 5 2 - 



AttLMb 



Vm 2 + q± 

^ L + C(e-^ L ), 



2M F 



0{e 



(A9) 



where the coupling constant g is defined as the vertex 
function T with all the particles being on the mass shell. 
In this way, the mass shift formula can be written as 



AM B (L) 



3g 2 



8ttLM b 



0{e 



-V2fiL 



)• (A10) 



This formula recovers Eq. (3.37) of Ref. [18| with to' = m 
and with symmetric factor 1/2 for identical particles. 



Note that to obtain positive 5 we need 



M B >m — m , 



(All) 



which is guaranteed by Eq. (|A1I) . All the above argument 
can be applied to the mass shift of to (to') through the 
m '-Mb loop (to- Mb loop), by replacing {Mb — > to, m — >■ 
Mb, to' — s- m'} ({Mb -> m' ,m -> Mb, m' -t to}). How- 
ever, in the small binding region (|AT|) . Eq. ([ATT]) is only 
valid for the self-energy of M B) so there is no pole con- 
tribution for the self-energies of intermediate particles of 
to and to'. The mass shift of the intermediate particles 
are then given by 



Ato =0{e' m ' L ) 
Am' =0(e- mL ) 

which do not alter the result (|A10|) . 



(A12) 
(A13) 
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